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I.  INTRODUCTION 


One  of  the  main  consequences  of  decreasing  the  yield  limit  in  a  threshold  test  ban  is  the 
increased  reliance  on  regional  distance  seismograms  for  monitoring  purposes.  Most  schemes 
for  determining  the  yield  at  regional  distances  are  based  on  measuring  the  amplitude  of  phases 
such  as  Pg  and  Lg.  It  is  extremely  important  to  a  TTBT  to  understand  the  nature  of  propa¬ 
gation  at  these  phases  in  different  regions,  as  well  as  possible  sources  of  bias  such  as  tec¬ 
tonic  release.  Here  we  report  on  our  continuing  research  on  tectonic  release. 

It  is  presently  thought  that  surface  waves  provide  the  best  data  for  determining  the 
orientation  of  the  equivalent  double  couple  representation  of  tectonic  release.  The  main  rea¬ 
son  for  this  is  that  until  recently  there  has  been  very  little  documentation  of  the  influence  of 
tectonic  release  on  P  waves.  The  most  widely  accepted  explanation  for  the  lack  of  an  obvious 
tectonic  release  signature  on  short-period  P  waves  is  that  the  tectonic  release  is  a  low  stress 
drop  phenomena  (Bache,  1976).  In  the  case  of  NTS,  the  low  stress  drop  will  conspire  with  the 
strike-slip  radiation  pattern  (which  does  not  radiate  P  waves  efficiently  to  teleseismic  distan¬ 
ces)  predicted  for  the  tectonic  release  to  make  the  short-period  signature  very  difficult  to 
observe  in  the  far-field  body  waves.  The  strongest  evidence  that  tectonic  release  can  radiate 
high  frequency  energy  is  presented  in  Wallace  et  al.  (1983,  1985),  where  it  is  shown  that  at 
regional  distances  the  long-period  P  waves  have  complexity  which  correlates  with  the  size  of 
the  SV  amplitude— suggesting  a  shear-dislocation  source.  Further,  the  SH  waves  can  be  mod¬ 
eled  with  a  strike-slip  orientation.  Recently  this  interpretation  of  the  P  wave  complexity  has 
been  reinterpreted  as  "slapdown”  (Douglas  et  al.,  1986).  In  Section  II  we  investigate  this 
possibility,  and  once  again  conclude  that  the  simplest  explanation  is  tectonic  release.  Further, 


the  long-period  tectonic  release  orientation  (surface  waves)  correlates  very  well  with  the 
azimuthal  variability  of  the  short-period  P  waves. 

There  is  an  apparent  frequency  dependence  of  the  tectonic  release  signature  at  regional 
distances.  At  periods  of  5-15  seconds,  Wallace  et  al.  (1983)  show  that  the  Pnj  waveforms  from 
Pahute  Mesa  explosions  are  sometimes  strongly  distorted  by  tectonic  release,  which  amounts  to 
superimposing  the  explosion  and  a  strike-slip  earthquake  source.  On  the  other  hand, 
Alexander  (1980),  Pomeroy  et  al.  (1982)  and  Gupta  and  Blandford  (1983)  present  results  which 
indicate  that  tectonic  release  has  little  effect  on  3-Hz  Lg  amplitudes.  A  similar  problem  is 
rectifying  the  far-field  representation  of  the  tectonic  release  with  the  observations  of  the 
tangential  component  of  strong  ground  motion.  The  observed  near-field  SH  waves  from  Pahute 
Mesa  events  are  much  more  complicated  and  significantly  smaller  than  expected.  Usually  the 
near-field  SH  waves  are  interpreted  in  terms  of  scattered  energy.  In  Section  III  we  discuss 
the  frequency  dependence  of  the  tectonic  release  signature  at  regional  distances.  The  only 
plausible  explanation  for  the  apparent  dispersion  is  a  shift  in  the  corner  frequency.  This  is 
consistent  with  the  results  discussed  in  Section  IV,  in  which  we  have  developed  a  "distributed- 
model  for  tectonic  release. 

Any  routine  inversion  procedure  developed  to  recover  tectonic  release  source  parameters 
will  require  the  analysis  of  sparse  data  sets,  requiring  the  use  of  waveforms  from  widely 
different  distances.  In  Section  V  we  present  a  procedure  for  the  joint  inversion  of  regional 
and  teleseismic  long-period  waveforms  for  source  parameter  estimates.  With  three  test  cases 
we  demonstrate  that  the  inclusion  of  regional  distance  seismograms  to  a  data  set  of 
teleseismic  observations  stabilizes  the  inversion  process  and  adds  valuable  constraints  to  the 
source  parameter  estimates,  particularly  when  the  teleseismic  Green’s  functions  are  inadequate. 
The  added  stability  and  constraints  from  regional  waveforms  can  be  attributed  to  the  robust 
nature  of  regional  distance  Green’s  functions,  as  well  as  the  added  coverage  of  the  focal 
sphere  that  regional  waveforms  provide.  Singular  value  decomposition  is  a  useful  inversion 


method  that  is  well  suited  for  this  study  because  it  allows  the  determination  of  the  influence 
that  individual  observations  are  having  on  the  model  parameter  estimates. 


II.  EVIDENCE  FOR  A  SHORT-PERIOD  TECTONIC  RELEASE  SIGNATURE 

In  two  earlier  papers  (Wallace  et  a!.,  1983,  1985),  we  discussed  the  evidence  for  tectonic 
release  from  underground  nuclear  explosions  on  Pahute  Mesa  at  the  Nevada  Test  Site  (NTS)  as 
observed  in  long-period  body  waves.  It  has  been  shown  for  some  time  that  the  nonisotropic 
component  of  the  surface  waves  from  most  of  these  events  could  be  explained  by  an 
equivalent  double-couple  source,  namely  strike-slip  motion  on  north-striking  faults.  In  an 
attempt  to  get  better  resolution  on  the  mechanism  of  the  tectonic  release,  we  modeled  the  SH 
wave  forms  using  the  constraints  imposed  by  reversals  in  the  SV  polarity  and  distortions  in 
the  long- period  P  waves.  We  found  that  the  orientation  inferred  from  the  body  waves  closely 
matched  those  determined  from  the  surface  waves.  Further,  we  were  able  to  determine  the 
seismic  moments  for  the  tectonic  release  by  comparison  with  several  shallow,  strike-slip 
earthquakes.  Although  the  SH  waveforms  provide  the  most  direct  measure  of  the  tectonic 
release,  the  tectonic  release  mechanism  adequately  explains  a  wide  range  of  long-period  body 
wave  observations.  In  particular,  we  identified  a  distortion  in  the  upper  mantle  long-period  P 
waves  which  we  attribute  to  the  phase  SP  from  the  tectonic  release.  The  amplitude  of  this 
phase  correlates  very  well  with  the  moment  of  the  tectonic  release;  it  can  be  modeled  by 
adding  a  strike-slip  synthetic  to  the  waveform  of  a  low  tectonic  release  event  (see  Figure  II- 
1);  and  finally  it  shows  a  change  in  polarity  (Figure  II-2)  as  expected  by  the  strike-slip 
radiation  pattern.  The  identification  of  the  sP  phase  provided  the  highest  frequency  signature 
of  tectonic  release  (0.33  Hz)  yet  known  and  suggested  that  short-period  P  waves  could  also  be 
affected  by  tectonic  release.  With  this  in  mind.  Lay  et  al.  (1984)  studied  the  first  cycle  (ab) 
amplitude  of  the  teleseismic  short-period  P  waves  from  the  Pahute  Mesa  events  and  found  a 
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Observations  at  LUB 


Figure  II-l :  A  comparison  of  the  P  and  PL  waves  for  BOXCAR  and  COLBY  at  LUB. 
Shown  below  is  the  COLBY  waveform  summed  with  a  synthetic  seismogram  to 
simulate  the  tectonic  release.  The  double  couple  has  a  pure  strike-slip 
orientation.  The  time  function  is  a  triangle  with  a  0.6  second  duration.  The 
seismic  moment  for  the  double  couple  required  to  obtain  the  fit  is  5  X  10^^ 
dyne-cm. 


Figure  I 1—2 :  A  comparison  of  the  GREELEY  waveforms  at  LUB  and  RCD.  Also  shown 
Is  the  waveform  of  the  8-16-66  earthquake.  This  earthquake  was  used  to 
simulate  the  GREELEY  record  by  summing  with  a  synthetic  explosion  waveform. 


systematic  azimuthal  variation  which  is  consistent  with  a  distortion  which  could  be  caused  by 
the  tectonic  release. 

There  is  some  controversy  over  our  interpretation  of  the  long-period  body  wave  tectonic 
release  signature.  Douglas  et  al.  (1986)  suggest  that  we  have  misinterpreted  the  complexity  in 
the  P  waves,  and  actually  the  distortion  is  caused  by  slapdown,  and  that  there  is  no  evidence 
to  support  that  tectonic  release  is  influencing  the  first  P  wave  arrivals  from  explosions  at 
NTS.  The  principal  data  Douglas  cites  are  the  broadband  recordings  at  EKA  of  the  large 
Pahute  Mesa  events  (Figure  II-3).  If  our  interpretation  of  tectonic  release  is  correct,  Douglas 
contends  that  the  deflection  he  has  marked  Ao  should  have  a  negative  polarity,  which  it  does 
not.  It  is  important  to  note  here  that  these  backswings  in  the  P  waves  are  interference 
patterns  between  the  explosion  P,  pP,  and  any  contribution  of  slapdown  with  that  of  the  sP 
phase  is  from  the  tectonic  release.  Such  interference  patterns  are  very  sensitive  to  the 
relative  timing  of  the  different  episodes.  A  small  shift  in  the  location  of  the  tectonic  release 
relative  to  the  explosion  can  produce  a  "misleading"  polarity.  A  direct  comparison  between  11- 
3  and  Figure  II -4  shows  that,  in  fact,  the  backswing  in  question  is  greatly  reduced  relative  to 
direct  P  for  the  high  tectonic  release  events  such  as  MUENSTER  and  BENHAM  at  5  EKA, 
which  is  consistent  with  the  tectonic  release  interpretation.  It  also  must  be  remembered  that 
SHA  (the  station  Wallace  et  al.  used  for  comparison)  is  24°  from  NTS,  while  EKA  is 
approximately  70°  away.  This  change  in  distance  corresponds  to  a  change  in  the  surface 
reflection  coefficient  of  sP  from  -0.55  to  -0.36;  or  a  reduction  by  one-third  at  EKA  as 
compared  to  SHA. 

We  want  to  reiterate  that  the  reason  we  favor  the  interpretation  that  tectonic  release 
can  have  a  signature  on  long-period  P  waves  is  its  consistency  with  a  large  data  set.  The 
amplitude  of  the  sP  phase  is  consistent  with  that  predicted  on  the  basis  of  the  SH  waves. 
The  SV  phase  from  BENHAM  and  GREELEY  shows  a  phase  reversal  which  is  consistent  with 
the  orientation  determined  for  the  tectonic  release  double  couple.  The  SV  reversal  is 


COLBY 


KASSERI 


FONTINA 


HANDLEY 


CAMEM8ERT 


GREELEY 


MUENSTER 


BENHAM 


WfW 


2  min. 


Figure  II-4:  The  vertical  component  records  for  8  large  Pahute  Mesa  explosions 
at  the  WWSSN  station  SHA.  These  are  long-period  seismograms  which  are  arranged 
according  to  the  importance  of  the  second  upswing  in  the  P  wave.  About  4j  min¬ 
utes  after  the  P  arrival  is  the  SV  wave. 


accompanied  by  a  reversal  in  polarity  of  the  phase  we  have  identified  as  sP.  The  long- period 
regional  distance  body  waves  (Pnl)  also  show  a  distortion  which  is  consistent  with  the 
interference  from  strike-slip  tectonic  release.  In  the  case  of  the  Pni  waveforms,  the  signature 
is  over  a  45-sec  window  rather  than  a  single  arrival  (sP).  Although  we  have  not  modeled 
slapdown,  and  as  Douglas  et  al.,  point  out,  it  can  have  an  effect  on  the  long-period  P  waves, 
we  believe  it  cannot  account  for  all  the  observations  listed  above.  In  particular,  slapdown 
cannot  produce  the  four-lobed  SH  radiation  pattern,  nor  the  polarity  reversal  of  sP.  Further, 
one  of  the  most  remarkable  facts  about  the  tectonic  release  from  NTS  is  its  consistency  in 
orientation.  With  only  two  or  three  exceptions  (HALFBREAK  and  PILEDRIVER),  the  tectonic 
release  is  predominantly  strike-slip  motion  on  north-trending  planes.  It  is  very  difficult  to 
appeal  to  a  form  of  explosion  source  anisotropy  or  slapdown  that  will  always  produce  an 
identical  apparent  double-couple  radiation  pattern. 

The  question  of  whether  the  tectonic  release  has  a  signature  at  the  frequencies  of  short- 
period  P  waves  is  a  separate  matter.  For  the  most  part,  the  short-period  waveforms  are  quite 
complicated,  and  no  one  has  convincingly  identified  an  equivalent  of  the  long-period  sP  phase. 
Lay  et  al.  (1984)  simply  plotted  the  corrected  ab  amplitude  of  25  Pahute  Mesa  explosions 
versus  azimuth  and  noted  that  they  appeared  to  have  a  sin  (2<f>)  variation.  Lay  et  al.  found 
that  the  best  fit  of  a  sin  {<(>)  curve  to  the  amplitudes  was  consistent  with  the  pattern  that 
would  be  expected  from  the  azimuthal  radiation  pattern  for  a  vertical  strike-slip  fault  with 
the  orientation  given  by  the  analysis  of  the  long-period  SH  waveforms.  The  significance  of 
the  sin  (4>)  fit  was  evaluated  with  an  f-test.  Nearly  all  the  Pahute  Mesa  explosions  pass  the 
test  at  the  99th  per  cent  level,  but  most  importantly,  the  events  with  the  largest  long-period 
tectonic  release  show  the  clearest  sin  (2^)  pattern.  There  are  alternative  explanations  for 
the  amplitude  behavior,  such  as  upper  mantle  structure  beneath  Pahute  Mesa  (Lynne  and  Lay, 
1984)  but  again  the  consistency  with  all  the  other  tectonic  release  data  suggest  that  the  ab 
amplitude  "could  be"  modulated  by  a  strike -slip  radiation  pattern. 


Lay  et  al.  also  cite  differences  in  waveform  between  stations  near  the  node  and  those  in 
the  lobe  as  evidence  of  tectonic  release.  Douglas  et  al.  suggest  that  the  change  in  waveform 
is  due  to  lateral  variations  in  attenuation.  In  particular,  they  map  the  loop  station  (UME) 
into  the  node  station  (OGD)  by  increasing  t*  by  a  factor  of  2  (see  Figure  4  in  Douglas  et  al.). 
One  way  to  test  this  path  bias  hypothesis  is  to  compare  UME  and  OGD  for  FAULTLESS, 
which  was  detonated  100  km  north  of  Pahute  Mesa;  the  waveforms  for  the  two  stations  are 
nearly  identical  (see  Lay  et  al.,  Figure  11).  Further,  we  can  test  the  stability  of  the  sin  (2</>) 
pattern  from  the  Pahute  Mesa  events  by  comparing  them  to  FAULTLESS.  Figure  7  in  Lay  et 
al.  shows  that  amplitude  pattern  from  FAULTLESS  is  significantly  different  from  that  of 
GREELEY;  in  fact,  they  are  nearly  reversed.  If  one  wants  to  appeal  to  attenuation  to  produce 
the  sin  (2 4)  pattern,  then  it  must  also  account  for  the  FAULTLESS  pattern.  It  seems  unlikely 
that  slight  travel  path  differences  between  FAULTLESS  and  GREELEY  to  teleseismic  distance 
stations  would  produce  the  factor  of  5  change  in  relative  amplitude. 

Recently,  Given  and  Mellman  (1985)  have  analyzed  the  surface  waves  from  several  of  the 
Pahute  events  using  sophisticated  phase  velocity  and  attenuation  correction  developed  for  NTS 
by  Stevens  et  al.  (1982).  Given  and  Mellman  determined  an  orientation  for  the  nonisotropic 
radiation  which  we  have  plotted  against  the  strike  of  sin  (2^)  pattern  observed  by  Lay  et  al. 
(1984)  for  a  set  of  eight  events  in  Figure  II-5.  The  strike  of  the  double  couple  varies  from 
90°  to  65°  for  the  different  events.  Although  there  is  some  difference  in  the  absolute  value 
of  the  strikes  as  determined  by  the  surface  and  body  waves,  their  agreement  is  good.  More 
remarkable  is  that  fact  that  as  the  strike  varies,  the  short  periods  track  the  long  periods  very 
well.  It  should  be  noted  that  this  correlation  breaks  down  for  events  outside  or  near  the 
edge  of  Silent  Valley  Caldera.  Again,  the  point  of  the  figure  is  that  the  short-period  P  waves 
amplitudes  correlate  with  the  long-period  tectonic  release.  Neither  structure  nor  attenuation 
can  be  ruled  out  as  causing  the  short-period  pattern,  but  there  is  no  reason  that  the  short- 
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Figure  II-5:  Strike  (in  degrees)  of  the  tectonic  release  double  couple  as 
determined  from  long-period  surface  waves  (triangular  symbols)  compared  with 
the  strike  of  the  sin  ( 2  4>  )  azimuthal  variation  observed  in  the  ab  amplitude  of 
the  short-period  P  waves.  All  of  the  events  are  within  silent  Canyon  Caldera. 
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period  P  waves  should  correlate  with  long-period  Love  and  Rayleigh  waves  unless  tectonic 
release  is  causing  the  modulation. 


One  difficulty  with  our  198S  paper  is  the  large  moments  we  determined  for  the  tectonic 
release.  These  values  were  determined  on  the  basis  of  comparison  of  SH  waves  from  two 
earthquakes:  (1)  16  August  1966,  near  the  Nevada-Utah  border,  and  (2)  15  March  1979,  in 
Homestead  Valley  in  south  central  California.  The  moment  values  were  very  consistent  from 
station  to  station  when  compared  to  either  earthquake.  Unfortunately,  the  average  values  for 
a  given  explosion  varied  by  a  factor  of  2  depending  on  the  reference  earthquake.  The  16 
August  1966  earthquake  gave  moments  which  were  twice  those  determined  using  the  15  March 
1979  event.  Since  we  did  not  know  which  estimate  was  more  reliable,  we  averaged  the  values 
to  produce  Table  I  in  the  1985  paper.  Since  that  publication,  several  investigators  have  shown 
that  there  is  a  strong  velocity  variation  across  the  Rocky  Mountain  front  (Grand  and 
Helmberger,  1984).  Since  the  16  August  1966  earthquake  is  along  the  front,  it  is  probably 
inappropriate  to  use  it  as  a  calibration  for  NTS.  Helmberger  et  al.  (1985)  have  developed  a 
lateral  varying  model  for  the  North  American  continent.  We  used  this  model  to  construct 
synthetic  seismograms  for  a  strike-slip  orientation  and  compare  them  to  GREELEY  in  Figure 
II-6.  As  can  be  seen,  the  waveform  fit  is  quite  good;  the  average  moment  we  obtain  is  1.9  X 
1024  dyne-cm,  which  is  in  very  close  agreement  with  the  value  we  determine  for  GREELEY  by 
the  comparison  with  the  15  March  1979  earthquake.  On  the  basis  of  this,  we  have  revised  the 
moments  for  tectonic  release  for  21  Pahute  Mesa  events.  The  new  values  are  given  in  Table 
II- 1. 

Since  our  revised  moments  are  smaller,  the  potential  for  tectonic  release  to  influence  the 
short-period  P  waves  is  also  reduced.  Although  we  are  confident  that  tectonic  release  can 
have  an  identifiable  signature  in  long-period  P  waves  (2  to  5  sec),  it  has  not  conclusively 
been  shown  that  it  influences  short-period  P  waves.  The  only  way  to  address  the  question  of 
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Figure  II-6:  A  comparison  of  the  SH  waveforms  from  the  explosion  GREELEY  with 
synthetics  generated  for  a  vertical  strike-slip  fault  with  the  orientation 
consistent  with  Pahute  Mesa  tectonic  release.  The  structure  model  used  is 
laterally  varying  (see  Helmberger  et  al.,  1985).  Above  the  observed  SH  pulse 
is  the  maximum  amplitude  in  10“3  cm,  while  the  number  below  the  synthetic  is 
the  moment  required  to  produce  the  observed  amplitude. 


TABLE  II- 1 


REVISED  MOMENTS  FOR  TECTONIC  RELEASE  FROM  PAHUTE  MESA  EXPLOSIONS 


Name 

Date 

(m-d-yr) 

Origin 

Time 

Lat. 

(°N) 

Long. 

(°W) 

Depth 

(km) 

mb 

i 

M0  X  1024 

Almendro 

06-06-73 

13:00 

37.24 

116.35 

1.064 

6.1 

0.61 

Benham 

12-19-66 

16:30 

37.23 

116.47 

1.402 

6.3 

3.43 

Boxcar 

04-26-68 

15:00 

37.29 

1 16.46 

1.158 

6.2 

0.86 

Camembert 

06-26-75 

12:30 

37.28 

116.37 

1.311 

6.1 

0.86 

Cheshire 

02-14-76 

11:30 

37.24 

116.47 

1.167 

5.8 

0.43 

Colby 

03-14-75 

12:30 

37.31 

116.47 

1.273 

6.2 

0.49 

Estuary 

03-09-76 

14:00 

37.31 

116.36 

0.869 

5.8 

o.43  ; 

Fontina 

02-12-76 

14:45 

37.27 

116.49 

1.219 

6.1 

1.35  1 

Greeley 

12-20-66 

15:30 

37.30 

116.41 

1.215 

6.3 

1.90 

Halfbeak 

06-30-66 

22:15 

37.32 

116.30 

0.819 

6.1 

0.61  ! 

Handley 

03-26-70 

19:00 

37.30 

116.53 

1.206 

6.4 

1.47  " 

Inlet 

11-20-75 

15:00 

37.22 

116.37 

0.817 

5.9 

0.18  j 

Jorum 

09-16-69 

14:30 

37.31 

1 16.46 

1.158 

6.1 

0.31 

Kusseri 

10-28-75 

14:30 

37.29 

116.41 

1.265 

6.2 

1.22  1 

Mast 

06-19-75 

13:00 

37.35 

116.32 

0.912 

5.9 

0.31 

Muenster 

01-03-76 

19:15 

37.30 

116.33 

1.451 

6.2 

1.96 

Pipkin 

10-08-69 

14:30 

37.26 

116.44 

0.617 

5.6 

0.12 

Pool 

03-17-76 

14:15 

37.26 

116.31 

0.879 

6.0 

0.12  | 

Scotch 

05-23-67 

14:00 

37.27 

116.37 

0.978 

5.7 

0.18  ' 

Stilton 

06-03-75 

14:20 

37.34 

116.52 

0.731 

5.8 

0.07  ‘ 

Tybo 

05-14-75 

14:00 

37.22 

116.47 

0.765 

5.9 

0.18  j 

III.  FREQUENCY  DEPENDENCE  OF  TECTONIC  RELEASE  AT  REGIONAL  DISTANCES 


At  regional  distances,  seismograms  are  very  complicated  due  to  the  waveguide  nature  of 
the  crust.  Wallace  (1986)  discusses  the  details  of  constructing  synthetic  seismograms  in  the 
distance  range  of  3°  to  14°.  The  long-period  wave  forms  for  explosive  sources  are  quite 
distinct  from  those  of  an  earthquake.  The  absence  of  S  in  the  source  and  the  high-frequency 
time  function  result  in  seismograms  that  do  not  have  well- developed  PL  phases.  On  the  other 
hand,  shallow  earthquakes  in  the  western  U.S.  with  magnitudes  greater  than  about  4.8 
routinely  write  measurable,  long-period  PL  phases  at  one  or  two  WWSSN  stations.  Under 
favorable  circumstances,  it  should  be  possible  to  separate  the  explosion  and  tectonic  release 
effects.  The  biggest  advantage  to  using  regional  distance  records  is  that  small  events  still 
produce  large  signal- to-noise  ratio  records  in  a  fairly  broad  frequency  band. 

Certain  systematic  affects  emerge  when  a  large  number  of  regional  distance  seismograms 
from  Pahute  Mesa  explosions  are  compared.  For  example,  the  Pm  waveforms  recorded  at  ALQ 
and  TUC  appear  much  more  like  those  produced  by  earthquakes  than  by  explosions.  Figure 
III- 1  is  a  composite  of  different  explosions  recorded  at  the  regional  WWSSN  network.  These 
long-period  seismograms  have  been  convolved  with  a  filter  whose  impulse  response  is  a 
triangle  with  a  2-sec  rise  and  fall.  Shown  below  the  observations  are  a  pair  of  synthetics 
which  have  been  similarly  filtered  (the  synthetics  were  generated  with  generalized  rays  and 
crustal  model  described  in  Wallace  et  al.,  1982).  These  synthetics  correspond  to  a  pure 
explosion  source  and  an  earthquake  source,  respectively. 

The  earthquake  source  has  a  strike-slip  mechanism  described  by  ToksOz  and  Kehrer 
(1972);  right-lateral  motion  on  a  plane  striking  N15°W.  ALG  and  TUC  are  in  the  radiation 
lobe  for  the  tectonic  release  PnJ,  while  LON,  which  is  in  good  agreement  with  the  explosion 


Figure  III-l:  The  Pn^  waveforms  for  several  Pahute  Mesa  explosions  (top  trace 
at  each  station)  and  synthetics  for  an  explosion  source  (middle  trace)  and  a 
double  couple  (bottom  trace).  Both  the  observations  and  synthetics  are  lightly 
filtered.  The  fault  orientation  for  the  double  couple  is  strike-slip;  the 
nodal  planes  are  sketched  on  the  location  map. 


synthetic,  is  near  a  node.  DUB,  which  is  in  the  negative  quadrant,  has  a  greatly  reduced  P 
amplitude.  Although  Figure  III- 1  is  a  composite  and  no  doubt  the  tectonic  release  varies  from 
event  to  event,  it  is  highly  suggestive  that  there  is  significant  tectonic  release  signature  on 
P„i  waveforms. 

FAULTLESS  was  detonated  about  100  km  north  of  Pahute  Mesa  at  Hot  Creek  Valley, 
Nevada.  In  contrast  to  the  Pahute  explosions,  the  FAULTLESS  Pnl  waveforms  show  much  less 
evidence  of  tectonic  release  associated  with  detonation  (FAULTLESS  does  appear  to  have 
numerous  strain  release  episodes  10  to  SO  seconds  after  detonation,  although  they  are  very 
small).  Figure  III-2  shows  the  regional  long-period  records.  Again  the  waveforms  have  been 
filtered  as  described  in  the  previous  paragraph.  Shown  below  each  observation  is  a  synthetic 
for  an  explosion  source.  Although  TUC  is  not  fully  explained  by  an  explosion  alone,  the 
tectonic  release  must  be  down  by  a  factor  of  3  in  comparison  with  GREELEY,  or  a  factor  of 
1.5  in  comparison  with  BOXCAR. 

In  an  attempt  to  isolate  the  tectonic  release  signature  on  Pn,  waveforms,  synthetics  were 
made  in  which  a  double  couple  was  added  to  an  explosion  synthetic.  The  attempt  here  is  to 
qualify  the  nature  of  tectonic  release  for  the  composite  records  shown  in  Figure  III- 1 ,  so  each 
record  is  fit  independently,  although  the  orientation  of  the  double  couple  was  constrained  to 
be  the  same  for  all  the  records.  The  explosion  synthetics  were  constructed  using  a  source 
time  function  described  by  Helmberger  and  Hadley  (1981).  They  used  their  source  time 
function  to  model  both  near-in  velocity  records  and  far-field  displacements  for  HANDLEY,  one 
of  the  explosions  considered  in  Figure  3;  so  the  values  they  determined  for  rise  time  and 
overshoot  (k  »  5,  B  =  2)  are  assumed  for  all  modeling.  Changing  the  values  of  k  and  B  within 
reasonable  limits  had  little  effect  on  the  filtered  synthetics.  This  is  similarly  true  if  we  had 
used  a  Haskell  (1967)  or  a  Von  Seggern  and  Bolandford  (1972)  source. 

Strike-slip  orientations  have  the  largest  effect  on  the  Pnl  waveforms  for  explosions.  The 
displacement  response  from  dip-slip  faulting  is  higher  frequency  than  that  for  strike-slip 
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Figure  III-2:  A  comparison  of  the  Pnj  waveforms  for  the  explosion  FAULTLESS 
(top  trace)  and  synthetics  computed  for  an  explosion  source.  Both  the 
observations  and  synthetics  have  been  lightly  filtered. 
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motion  (this  is  a  result  of  excitation)  and  generally  adds  or  subtracts  to  the  explosion 
waveform  without  substantially  changing  it.  Since  the  ration  of  dip-slip  to  strike-slip  motion 
is  small  on  the  basis  of  Wallace  et  al.  (1984),  only  the  orientation  of  the  strike-slip  component 
was  determined.  Different  orientations  for  the  strike  were  tested  for  compatibility  with  the 
observations.  Figure  III-3  summarizes  the  analysis.  Shown  are  the  same  explosions  as  in  ■ 

i 

Figure  III-l,  and  synthetics  generated  for  a  combination  of  double  couple  and  explosion.  j 

The  LON  and  COR  records  have  a  profound  effect  on  the  strike  of  the  double  couple.  \ 

( 

The  contribution  of  the  explosion  to  the  waveform  is  much  greater  than  that  of  the  double 
couple,  implying  that  these  stations  are  near  the  node  of  the  radiation  pattern.  This  is 
particularly  true  of  LON  for  nine  different  explosions.  The  best  fitting  strike-slip  fault  has  a  < 

strike  of  N20°W  (not  significantly  different  from  that  of  ToksOz  and  Kehrer  (1972).  The  ■ 

numbers  to  the  right  of  each  seismogram  pair  in  Figure  III-3  give  the  ratio  of  Pn  displacement  j 

caused  by  the  explosion  contribution  to  that  caused  by  the  strike-slip  dislocation.  Since  ; 

I 

different  explosions  are  used  in  the  composite  figure,  no  absolute  moments  are  given.  Rather,  ' 

l 

the  ratios  give  a  measure  of  the  importance  of  tectonic  release  for  a  given  azimuth.  J 

Figure  III-4  shows  the  vertical  component  for  12  Pahute  explosion  waveforms  at  PAS,  a  j 

long-period  (30-90)  station  5.5°  from  NTS.  At  the  periods  of  these  records,  it  is  assumed  that  ^ 

the  travel  paths  are  identical,  and  most  of  the  differences  in  the  waveform  must  be  due  to 
source.  It  is  interesting  to  compare  events  which  are  spatially  very  close  but  separated  in 
time.  For  example,  GREELEY  (12-20-66)  and  KASSERI  (10-28-75)  are  less  than  3  km  apart 
and  about  the  same  yield,  but  GREELEY  has  a  much  higher  tectonic  release.  In  Figure  III-4  | 

the  PAS  records  are  very  similar,  except  that  KASSERI  is  deficient  in  the  long  periods.  j 

Single  station  comparisons  of  many  events  should  make  it  possible  to  develop  an  algorithm  to 
identify  source  effects  such  as  depth  of  burial  as  well  as  tectonic  release. 

It  is  a  fairly  straightforward  process  to  invert  for  the  moment  tensor  (MT)  of  a  seismic 
source,  provided  adequate  Green’s  functions  can  be  constructed,  although  for  shallow 
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Figure  III-3:  A  comparison  of  the  Pn^  waveforms  of  the  same  explosions  as  in 
Figure  I I I— 1  and  synthetic  explosion  waveforms  with  a  component  of  tectonic 
release.  Both  observations  and  synthetics  have  been  lightly  filtered.  The 
numbers  to  the  right  of  each  seismogram  pair  give  the  ratio  of  Pn  amplitude  due 
to  the  double  couple  and  that  due  to  che  explosion. 
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Figure  III-4:  Vertical  seismograms  for  12  Pahute  Mesa  explosions  at  PAS  (30- 
90).  The  station  is  3.3°  from  NTS.  Numbers  to  right  of  traces  give  the 
maximum  amplitude  in  mm. 


. 


explosions',  high-frequency  data  are  required  to  constrain  the  Mci  components.  The  explosion 
moment  is  the  average  of  the  trace  of  the  moment  tensor,  and  the  deviatoric  component  is 
the  remainder.  We  have  inverted  the  regional  body  waves  recorded  on  the  broadband 
Lawrence  Livermore  National  Laboratory  network  for  three  Pahute  Mesa  explosions  (CHES1RE, 
FARM  and  POOL)  in  different  frequency  bands  in  an  attempt  to  isolate  any  frequency 
dependence  of  the  deviatoric  component  of  the  MT. 

At  long  periods  (>2  seconds)  a  two-layered  crust  over  a  halfspace  mantle  was  a  sufficient 
travel  path  structure  to  recover  the  moment  tensor.  Table  III- 1  shows  the  isotropic  and 
major  double  couple  components  of  the  MT.  Also  shown  is  the  size  of  the  minor  double 
couple,  an  indication  of  the  likelihood  that  the  major  double  couple  actually  represents 
tectonic  release.  The  inversion  was  performed  in  a  linear  fashion,  which  requires  that  all  the 
elements  of  the  MT  have  the  same  time  history.  It  was  assumed  that  the  explosion  had  a 
Helmberger-Hadley  time  history  with  k  =  15,  B  =  0.  Although  this  is  a  poor  representation 
for  the  tectonic  release  time  history,  it  is  probably  sufficient  considering  the  uncertainty  in 
the  travel  path  structure. 

Figure  III-5  shows  the  frequency  dependence  of  the  deviatoric  part  of  the  moment  tensor 
for  the  three  NTS  events.  The  LLNL  records  were  band-pass  filtered  (the  high-frequency  cut¬ 
off  is  shown  on  the  abscissa).  The  resulting  moment  tensor  had  the  long-period  isotropic 
component  subtracted  from  it,  and  the  largest  remaining  eigenvalue  was  assumed  to  represent 
the  tectonic  release  (the  direction  of  stress  axes  was  not  considered).  As  can  be  seen,  the 
apparent  size  of  the  tectonic  release  decreases  with  increasing  frequency.  It  is  difficult  to 
explain  the  frequency  dependence  of  the  tectonic  release.  At  long  periods,  the  regional- 
distance  data  is  consistent  with  the  teleseismic  analysis. 

Murphy  et  al.  (1986)  studied  the  mb  anomaly  associated  with  the  RULISON  explosion  (10 
September  1969,  yield  40  kt).  The  observed  mb  is  up  to  0.3  units  smaller  than  expected  based 
on  a  comparison  with  GASBUGGY.  Murphy  et  al.  attribute  this  anomaly  to  the  signature  of 
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Figure  III-5:  The  ratio  of  the  deviator ic  to  explosion  moment  for  Pahute  Mesa 
explosions  CHESHIRE,  POOL  and  FARM.  M  is  only  determined  at  long  periods. 


tectonic  release.  In  this  case,  the  tectonic  release  orientation  is  a  normal  fault,  so  it 


subtracts  from  the  compressional  P  wave  arrivals.  Murphy  et  al.  attribute  this  anomaly  to  the 


signature  of  tectonic  release.  In  this  case,  the  tectonic  release  orientation  is  a  normal  fault. 


so  it  subtracts  from  the  compressional  P  wave  arrivals.  Murphy  et  al.  further  state  that  the 


amplitude  of  the  vertical  component  Lg  is  affected  by  the  tectonic  release,  thus  still 


providing  a  good  yield  estimator. 


i 


VI.  NEAR-FIELD  TECTONIC  RELEASE 


The  tangential  motions  recorded  in  the  near-field  from  underground  nuclear  explosions 
are  inconsistent  with  the  teleseismic  representation  of  the  tectonic  release.  As  an  example  of 
this,  Figure  VI- 1  shows  the  three  components  of  ground  motion  from  BOXCAR  recorded  at  a 
distance  of  7.4  km.  Also  shown  in  Figure  VI- 1  are  synthetics  for  an  explosion  source  function 
and  a  double  couple.  The  velocity  structure  that  was  used  for  the  synthetics  is  that 
determined  by  Barker  et  al.  (1985).  The  P  wave  arrivals  are  matched  to  first  order  by  the 
explosion  synthetic.  The  SH  waves  were  generated  using  a  moment  of  0.86  X  1024  dyne-cm 
(Wallace  et  al.,  1986).  As  can  be  seen,  the  tangential  component  is  poorly  matched.  This 
gross  mismatch  indicates  two  things:  (1)  significant  off-azimuth  energy,  and  (2)  the 
inadequacy  of  a  point  source  representation  for  the  tectonic  release.  On  the  basis  of  our 
teleseismic  modeling  of  the  SH  energy  from  nuclear  explosions,  we  would  argue  strongly  that 
there  is  appreciable  release  of  tectonic  strain,  and  that  the  tangential  energy  is  not  a  product 
of  asymmetry  of  the  explosion  source  (see  Wallace  et  al.,  1985,  1986).  If  the  tectonic  release 
is  purely  a  triggered  fault  motion  as  suggested  by  Aki  and  Tsaiu  (1972),  then  it  would  be 
expected  that  the  strong-motion  SH  waves  would  be  similar  to  those  produced  by  a  similar¬ 
sized  earthquake.  Figure  VI-2  compares  the  SH  waves  from  an  aftershock  of  the  1979 
Imperial  Valley  earthquake  (Ml  »  5.0,  Mo  **  0.7  X  1024  dyne-cm)  recorded  at  10.4  km  distance 
with  the  BOXCAR  record  and  a  synthetic  for  a  simple  dislocation  model.  Apparently,  the 
point  source  representation  is  sufficient  for  the  earthquake,  implying  that  at  least  some 
component  of  the  tangential  record  from  the  explosion  is  due  to  a  mechanism  other  than 
triggered  fault  motion. 
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Figure  IV-1 :  Comparison  of  strong  ground  motion  from  BOXCAR  observed  at  a 

recording  site  7.4  km  away  from  the  source  with  synthetic  seismograms 
calculated  using  the  model  of  Hartzell  et  al.  (1983).  The  tangential  model  is 
from  the  best  fitting  point-source  double  couple. 
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Fij’ure  I V— 3 :  Distributed  source  SH  waves  in  the  near  and  far  field.  Tin* 
scroll)'  motion  fit  requires  a  total  moment  of  0.18  X  10^'  dyno-cm.  The  la r- 
fleld  comparison  is  for  a  single  point  source  with  a  moment  of  0.18  X  10^'  and 
16  sources  of  0.22  X  10^*  dyne -cm.  There  Is  no  difference  in  waveform  or 
amp  I i tude . 
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Kl««re  IV-2:  Comparison  of  BOXCAK  and  an  Imperial  Valley  af tersluiek  with  a 
svntlietlc  SH  wavetraln.  I  he  explosion  and  earthquake  have  approximately  the 


ar-lleld  moment. 
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Archambeau  (1972)  proposed  a  mechanism  for  tectonic  release  which  is  due  to  the  loss  of 
strength  of  the  material  in  a  volume  about  the  explosion.  In  an  attempt  to  simulate  this 
model,  we  constructed  a  spherical,  symmetric,  distributed  source.  We  placed  16  point  sources 
(of  equal  moment)  about  the  working  point.  The  elastic  radius,  or  extent  of  shear  stress 
drop,  expected  for  Archambeau’s  model  is  given  by: 

Re  *  FpasV>{oo)  Act 

where  F  is  the  ratio  of  the  tectonic  release  moment  to  the  explosion  moment.  Act  is  the  stress 
drop,  p  and  a  are  the  source  region  density  and  P  velocity,  respectively,  and  \l>(oo)  is  the  static 
level  of  the  explosion  displacement  potential.  For  BOXCAR  (1000  kt)  and  the  F  factor 
calculated  on  the  basis  of  the  tectonic  release  moment  of  0.87  X  1024  dyne-cm,  the  elastic 
radius  should  be  on  the  order  of  1.1  km  if  the  stress  drop  is  100  bars.  This  radius  seems 
extraordinarily  large,  and  indeed  would  have  a  significant  effect  on  a  recording  only  7  km 
away.  The  value  of  100  bars  for  stress  drop  is  taken  from  the  "average"  intraplate  earthquake 
(Liu  and  Kanamori,  1983).  For  earlier  work  (Wallace  et  al.,  1983),  it  was  shown  that  the  sP 
waveforms  from  BOXCAR  required  a  much  higher  stress  drop:  30-500  bars.  This  is  consistent 
with  the  work  of  McKeown  and  Dickey  (1969),  who  showed  that  the  aftershock  distribution  of 
explosions  with  large  F  factors  is  confined  to  a  region  much  smaller  than  expected  for  an 
earthquake  of  comparable  moment.  Using  a  stress  drop  of  500  bars,  the  elastic  radius  for 
BOXCAR  is  0.3  km.  We  used  this  radius  to  define  the  sphere  of  the  distributed  source. 

Figure  IV -3  shows  a  comparison  of  the  observed  tangential  and  distributed  source 
synthetic  waveforms  for  BOXCAR.  Although  the  match  is  not  perfect,  it  does  preserve  the 
character  of  the  interference  pattern.  Also  shown  in  Figure  IV-3  is  the  long-period 
teleseismic  SH  pulse  from  the  distributed  source  compared  with  that  from  a  point  source.  It 
is  apparent  that  the  distributed  source  model  will  satisfy  the  local  data  and  far-field 
observation,  although  the  seismic  moment  that  is  required  in  the  near  field  is  one  half  that 
derived  on  the  basis  of  teleseismic  data  alone. 
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One  of  the  main  problems  in  dealing  with  the  near-field  data  is  the  poor  azimuthal  and 
range  coverage.  It  is  difficult  to  see  coherence  in  the  SH  pulse  or  node  changes.  One 
exception  appears  to  be  HALFBEAK  (30  June  1966,  yield  300  kt).  This  event  has  a  large  F 
factor  (1.9,  using  the  moment  determined  by  Wallace  et  al.,  1986)  and  has  very  coherent  SH 
arrivals  in  the  near  field.  Figure  I V-4  shows  the  tangential  components  arranged  in  distance 
from  0.91  to  3.0  km.  Stations  4,  5,  and  6  are  in  a  line  striking  north  of  HALFBEAK.  This 
line  of  stations  is  along  a  direction  predicted  to  be  a  maximum  for  the  tectonic  release  SH 
(on  the  basis  of  the  teleseismic  modeling).  Station  S8  is  to  the  east,  and  is  predicted  to  have 
a  polarity  reversal  as  compared  to  S6.  On  the  basis  of  timing,  this  reversal  can  be  seen, 
although  it  is  not  as  dramatic  as  expected. 

Figure  IV -5  shows  the  tangential  recordings  from  HALFBEAK  and  synthetics  generated 
for  a  distributed  shear  dislocation  source.  The  velocity  structure  that  was  used  for  the 
synthetics  is  that  determined  by  Barker  et  al.  (1985).  Green’s  functions  were  computed  with 
Kennet-Bouchen  numerical  integration.  The  source  orientation  for  the  shear  dislocations  was 
determined  by  a  constrained  moment  tensor  inversion.  It  was  assumed  that  the  orientation 
was  strike-slip,  and  the  tangential  waveforms  were  inverted  for  strike.  The  resulting  value, 

8  =  N7°E,  is  consistent  with  larger  Pahute  Mesa  events  studied  teleseismically  (Wallace  et  al., 
1985,  1986).  The  inversion  produced  a  stable  estimate  of  strike,  but  the  quality  of  fit  to  the 
waveforms  was  quite  poor.  Since  the  near-field  recordings  are  only  2-4  times  the  depth  of 
burial  (d  «  815  m),  tectonic  release  from  a  distributed  source  will  produce  a  very  complicated 
interference  pattern.  In  an  attempt  to  simulate  this  effect,  we  placed  16  point  sources  (of 
equal  moment,  the  total  of  which  was  determined  by  the  moment  tensor  inversion)  about  the 
working  point  at  the  distance  of  1  elastic  radius  (given  by  Archambeau,  1972,  and  a  stress 
drop  of  100  bars).  Although  the  fit  in  Figure  IV-5  is  not  perfect,  it  does  offer  a  reasonable 
explanation  for  the  data.  In  the  next  step,  we  plan  to  use  this  model  for  tectonic  release  to 
fit  the  vertical  and  radial  near-field  observations. 
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Figure  IV-A:  Tangential  Displacements  for  HALFBEAK.  Dashed  outline  is  the  Si¬ 
lent  Canyon  Caldera.  The  S4-S6  profile  is  near  the  SH  max. 
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Figure  IV-5:  A  comparison  of  the  observed  and  synthetic  near-field  records 
from  HALFBEAK.  .  Top  trace  is  the  observed  velocity.  The  synthetics  were 
calculated  with  a  distributed  shear  dislocation  source. 


V.  JOINT  INVERSION  OF  REGIONAL  AND  TELESEISMIC  WAVEFORMS 


The  inconsistency  of  data  from  different  distances  makes  the  absolute  determination  of 
tectonic  release  source  parameters  suspect.  Any  tectonic  release  model  must  explain  the 
discrepancy  between  regional  and  teleseismic  estimates  of  the  tectonic  release  size. 
Unfortunately,  many  explosions  contain  only  a  few  usable  teleseismic  records,  so  in  general  we 
are  dealing  with  a  sparse  data  set.  We  have  developed  a  procedure  for  the  joint  inversion  of 
regional  and  teleseismic  long-period  waveforms  for  source  parameter  estimates.  At  present, 
the  procedure  has  only  been  used  on  earthquake  waveforms. 

The  source  parameters  of  moderate  to  large-sized  earthquakes  are  widely  used  to  inter 
the  present  day  tectonic  environment.  In  most  regions,  large  earthquakes  are  infrequent,  and 
moderate-sized  events  must  be  used  to  investigate  the  style  of  deformation  and  state  of  stress. 
Unfortunately,  useful  recordings  of  moderate-sized  events  at  teleseismic  distances  are  often 
sparse,  making  it  difficult  to  constrain  the  source  mechanism  through  the  inversion  of 
teleseismic  long-period  body  waves  alone.  Wallace  et  al.  (1981)  introduced  a  method  for  the 
determination  of  source  parameters  by  the  inversion  of  regional-distance  long-period  body 
waves.  This  method  has  proved  valuable,  and  as  little  as  three  stations  (azimuthally 
distributed)  are  required  in  order  to  constrain  a  solution.  The  reason  for  this  is  the  nature 
of  Pn[  wave  trains:  the  combination  of  P  and  SV  rays  provides  a  greater  coverage  of  the 
outer  portion  of  the  focal  sphere  than  teleseismic  P  and  SH  waves,  which  have  much  steeper 
take-off  angles.  In  addition,  as  much  as  90  seconds  of  record  can  be  fit  in  regional  waveform 
inversions,  whereas  most  useful  information  from  shallow-source  mechanisms  occurs  in  the  first 
30  seconds  of  the  teleseismic  waveform.  On  the  other  hand,  regional  waveforms  are  relatively 
insensitive  to  source  depth,  and  teleseismic  information  can  be  used  to  constrain  this 


parameter.  Ideally,  a  method  which  inverts  both  of  these  data  simultaneously  will  provide  the 
best  constraints  on  source  parameters.  A  joint  inversion  of  teleseismic  P  and  SH,  and 
regional  vertical  and  radial  component  long-period  waveforms,  has  the  distinct  advantage  of 
examining  energy  with  take-off  angles  that  range  from  18°-52°,  considerably  better  coverage 
than  teleseismic  waveform  or  regional  waveform  inversion  alone.  This  paper  presents  a 
technique  for  the  simultaneous  inversion  of  teleseismic  and  regional  long-period  waveforms  for 
source  parameters.  The  utility  and  benefits  of  the  inversion  method  will  be  demonstrated  on 
synthetically  derived  data  for  three  cases. 

The  Forward  Model 

The  teleseismic  and  regional  P  and  S  wave  synthetics  can  be  expressed  in  cylindrical 
coordinates  as 

s 

Wp(r,z,0,t)  =  Mo/4*p0(S(t)*  £  G,(t)At) 

i=i 

2 

Wx(r,z,0,t)  -  Mo/4*p0(S(t)*  £  Gj(t)Aj+3) 

j=t 

where  wp  and  Ws  are  the  P  wave  and  SH  wave  displacements,  respectively,  Mc  is  the  seismic 
moment,  p0  is  the  source  region  density,  and  S(t)  is  the  far-field  time  history.  The  Afs  are 
the  source  coefficients  determined  by  source  orientation: 

Aj(0,A ,6)  *  sin  26  cos  A  sin  6  +  (l/2)cos  26  sin  A  sin  25 
A2(6,A,6)  =  cos  6  cos  A  cos  6  -  sin  6  sin  A  cos  26 
A3(6,X,6)  =  (l/2)sin  A  sin  26 

where  6  is  the  receiver  azimuth  from  the  end  of  the  fault  plane,  A  is  the  rake  angle,  and  6  is 
the  dip  angle. 

The  G;  terms  represent  the  Green’s  functions  for  the  three  fundamental  faults;  = 


vertical  strike-slip,  G2  =  vertical  dip-slip,  and  G3  =  45  dip-slip.  The  Green’s  functions  for  the 
regional  distance  are  for  either  the  vertical  or  radial  component  and  were  generated  using 
generalized  ray  theory. 
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We  employed  the  inversion  method  described  by  Wallace  et  al.  (1981),  which  uses  the 
error  function 

e  =  1  -  Jfg/(Jf2)1/2(Jg2)1/2  (1), 

where  f  is  the  observed  record  and  g  is  the  synthetic  waveform.  There  is  an  error  associated 
with  each  waveform  observation  and  the  sum  of  the  squares  of  these  errors  can  be  minimized 
with  respect  to  the  three  fault  parameters  strike,  dip,  and  rake  (  9,  S  and  A)  through  an 
iterative  nonlinear  inversion.  The  generalized  inverse  for  this  inversion  can  be  expressed  in 
the  form  of  singular  value  decomposition  (SVD)  as  VpA"£Up,  where  Vp  is  a  p  X  p  matrix  of 
model  eigenvectors  Vj  (p  =  number  of  nonzero  single  values)  and  Up  is  a  N  X  p  matrix  of  data 
eigenvectors  U;  (N  =  number  of  observations).  The  term  A'p  is  a  p  X  p  matrix  composed  of 
diagonal  elements  of  the  form  1/A;,  where  A;  are  the  eigenvalues  or  single  values.  The  model 
parameter  vector  change  can  thus  be  expressed  as  the  weighted  sum  of  the  model  eigenvectors 

p 

as  £  ajVi,  where 
i=l 

a*  «  l/AiUiTe  (2). 

One  benefit  of  SVD  is  that  the  influence  a  single  station  or  set  of  stations  are  having  on  the 
model  parameter  estimates  can  be  easily  determined.  To  increase  stability  in  the  nonlinear 
inversion,  a  damping  factor  can  be  included  so  that  the  generalized  inverse  is  in  the  form  Vp 
diag(A;/Af  +  0)Uj,  where  0  is  the  damping  factor.  This  takes  care  of  the  problem  of  small 
single  values,  which  can  cause  unusually  large  changes  in  the  model  parameter  estimates 
during  the  iteration  process.  The  maximum  likelihood  inverse  or  weighted  least-squares  can 
also  be  used  when  it  can  be  determined  that  noisy  observations  are  imposing  a  strong 
influence  on  the  model  parameter  estimates. 

Results 

Three  test  cases  were  studied  to  illustrate  the  strengths  of  joint  inversion  of  regional 
and  teleseismic  data.  The  first  case  is  a  pure  strike-slip  mechanism  shown  in  Figure  V-l. 
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There  were  eight  teleseismic  observations  (four  P  wave  and  four  SH)  and  four  regional 
stations  with  vertical  and  radial  components  from  each.  The  synthetic  data  were  derived 
using  a  source  depth  of  10  km.  Both  the  synthetic  data  and  the  Green’s  functions  for  the 
teleseismic  distances  were  generated  using  the  three  rays,  P,  pP  and  sP.  The  starting  solution 
was  a  pure  normal  fault  (6  =  130°,  6  =  45°,  and  A  =  -90°).  Figure  V-l  shows  the  solution  at 
3,  6  and  10  iterations.  The  exact  solution  was  obtained  after  12  iterations.  Also  shown  in 
Figure  V-l  is  a  plot  of  the  error  function  values  from  each  inversion  versus  the  Green’s 
function  source  depth.  A  minimum  error  of  zero  was  obtained  using  the  10-km  Green’s 
functions,  thus  constraining  the  source  depth  exactly. 

Table  V-l  contains  the  results  for  the  same  problem  except  that  an  inappropriate  time 
function  was  used.  Table  I  lists  the  solutions  from  each  inversion  using  Green’s  functions 
generated  for  depths  between  5-16  km.  These  results  indicate  that  the  inversion  for  model 
parameter  estimates  is  rouust  for  Green’s  function  depths  of  6-12  km,  with  strike,  dip  and 
rake  only  deviating  from  the  true  solution  by  l°-3°.  The  error  functions  were  minimized  at  7 
km  (true  depth  =  10  km).  Therefore,  the  time  function  does  not  strongly  influence  the  model 
parameter  estimates  strike,  dip  and  rake  in  the  simultaneous  inversion;  however,  it  does 
influence  the  depth  estimate.  In  this  example,  the  long  time  function  produced  an 
underestimation  of  the  source  depth. 

The  next  case  addresses  the  problem  of  inadequate  Green’s  functions.  This  case  has  the 
same  station  distribution  as  the  previous  example.  The  teleseismic  synthetic  data  were 
generated  using  13  rays,  which  consisted  of  many  combinations  of  crustal  bounces.  The 
effects  of  a  complicated  receiver  structure,  which  had  50  rays  with  many  combinations  of 
crustal  reverberations  and  P-SV  conversions  in  a  5-layer  crust,  were  also  included.  The 
effects  of  P-SV  conversions  on  the  vertical  component  are  small  due  to  the  steep  ray  path  for 
teleseismic  P-wave  energy. 
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Figure  V-l :  Source  mechanism  and  waveform  fits  at  3,  6,  and  10  iterations. 
The  error  (equation  1)  is  next  to  each  waveform  (not  all  waveforms  used  in  the 
inversion  are  shown).  The  synthetic  observed  waveform  is  on  top,  and  the 
synthetic  predicted  waveform  is  on  the  bottom.  The  error  for  solutions 
obtained  using  different  source  depths  is  also  shown.  True  source  depth  was 
10  km. 
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TABLE  V-l 


INVERSION  RESULTS  USING  SUITE  OF  GREEN’S  FUNCTIONS 
AND  AN  INAPPROPRIATE  TIME  FUNCTION4 


Depth  (km) 

Errorb 

Strike 

Dip 

Rake 

5 

.2344 

211.6 

88.5 

2.5 

6 

.2037 

34.6 

89.6 

-2.0 

7 

.1271 

30.7 

88.4 

-0.7 

9 

.2159 

32.5 

88.9 

-3.2 

10 

.2093 

32.6 

86.7 

-0.6 

11 

.3374 

30.7 

85.6 

-3.0 

12 

.2410 

30.8 

86.4 

2.3 

14 

.4914 

30.7 

75.8 

10.9 

16 

.4745 

41.8 

71.0 

3.6 

“Observed  synthetic 

data  had 

a  trapezoid 

time  function 

with  rise  time. 

rupture  time,  and  fall  time  =  1,  2,  1. 


bAverage  size  of  the  error  for  all  observations.  Error  is  defined  in  equation 
1  in  text. 
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The  results  of  the  inversion  using  teleseismic  data  alone  for  this  oblique-slip  case  are 
shown  in  Figure  V-2a.  Again,  the  teleseismic  Green’s  functions  were  generated  using  only  3 
rays  and,  therefore,  are  inadequate  to  completely  model  the  observations.  The  true  depth  was 
10  km,  but  the  error  functions  were  minimized  using  Green’s  functions  generated  for  a  6-km 
depth.  Table  V-2  summarizes  the  source  parameters  obtained  from  the  inversion  using  the 
suite  of  Green’s  functions.  The  error  forms  a  broad  minimum  between  6  and  9  km. 
Therefore,  the  source  depth  from  this  inversion  is  poorly  constrained.  The  solution  from  the 
teleseismic  data  alone  (6  =  60°,  6  *  50°,  A  =  66°)  is  significantly  different  from  the  true 
solution,  yet  the  waveform  fits  look  very  reasonable. 

Figure  V-2b  shows  the  inversion  results  when  the  regional  waveforms  are  included.  The 
regional  distance  synthetic  data  for  this  example  were  generated  using  a  crustal  thickness  of 
30  km.  The  Green’s  functions  were  generated  using  a  crustal  thickness  of  40  km  and  are, 
therefore,  also  inadequate.  The  inclusion  of  the  regional  waveforms  in  the  inversion,  however, 
had  a  significant  effect  on  constraining  the  predicted  solution  to  the  true  solution,  even 
though  the  regional  distance  Green’s  functions  are  inadequate.  Table  V-3  shows  the  results  of 
including  the  regional  waveforms  in  the  inversion.  For  all  depths,  although  the  size  of  the 
error  might  change,  the  model  parameter  estimates  are  within  l°-5°  of  the  true  solution.  The 
addition  of  regional  waveforms  to  the  data  set,  therefore,  has  a  significant  stabilizing  effect 
on  the  inversion.  This  result  demonstrates  the  robustness  of  the  regional  distance  Green’s 
functions  and  is  consistent  with  the  findings  of  Wallace  (1986)  that  regional  distance 
waveforms  are  relatively  insensitive  to  the  fine  details  of  crustal  structure  and  crustal 
thickness,  thus  making  them  ideal  for  source  parameter  studies. 

The  teleseismic  distance  station  distribution  (seven  P  waves  and  four  SH)  and  the  pure 
normal  mechanism  of  the  third  case  are  shown  in  Figure  V-3a.  This  example  suffers  from 
poor  teleseismic  azimuthal  coverage,  a  typical  problem  encountered  for  moderate-sized  events. 
With  the  teleseismic  data  alone,  the  inversion  suffered  from  nonuniqueness  and  instability  from 
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Figure  V-2a:  Inversion  results  from  teleselsmic  data  alone  using  inadequate 
teleseisraic  Green's  functions.  Note  the  acceptable  waveform  fits  with  a 
solution  that  deviates  significantly  from  the  true  mechanism. 

Figure  V-2b:  Results  from  joint  inversion  of  teleseismic  and  regional 
waveforms.  Regional  waveforms  were  added  to  the  teleseismic  data  set  from 
Figure  V-2a.  Both  teleseismic  and  regional  distance  Green's  functions  were 
inadequate.  The  Inclusion  of  the  regional  data  in  the  inversion  produced  a 
solution  very  close  to  the  true  mechanism,  even  though  the  regional  distance 
Green's  functions  are  inadequate.  Stations  1  and  3  are  near  a  nodal. 
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INVERSION  RESULTS,  TELESEISMIC  DATA  ALONE, 
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Figure  V-3a:  Solution  space  for  pure  normal  fault  mechanism  with  poor  station 
distribution.  The  average  size  of  the  error  function  is  contoured  over  the 
variation  of  strike  and  rake  while  dip  is  held  constant  at  70°.  Elongated  flat 
low-misfit  region  around  true  solution  indicates  a  poorly  resolved  mechanism. 


Figure  V-3b:  Solution  space  for  normal  fault  after  adding  four  regional 

stations  (boxes).  Mechanism  is  much  better  constrained  and  elongated  flat 
region  has  disappeared. 


small  eigenvalues.  Using  a  damping  factor  of  0.1  and  a  starting  solution  close  to  the  true 
solution  (6  *  40°,  5  =  70°,  A  =  -90°),  the  inversion  arrived  at  a  nonunique  solution  (9  =  -48°, 
S  =  25°,  A  =  -52°).  Using  a  much  larger  constant  damping  factor  of  3.0,  the  solution  did  not 
deviate  far  from  the  true  solution;  however,  after  75  iterations,  the  parameter  estimates 

(9  a  65°,  S  *  61°,  A  =  -81°)  had  still  failed  to  converge  to  the  correct  solution.  The  reason 

for  this  behavior  can  be  seen  in  Figure  V-3a.  This  figure  shows  the  solution  space  for  this 
problem.  The  average  length  of  the  error  function  is  contoured  over  the  variation  of  strike 
and  rake  while  dip  is  held  constant  at  70°.  An  elongated  trough  forms  a  region  of  low  misfit. 
The  true  solution  lies  close  to  the  center  of  this  trough,  which  explains  the  slow 
convergence  and  indicates  nonuniqueness  (anywhere  in  this  region  the  solution  is  satisfactory). 

With  the  addition  of  data  from  the  four  regional  stations,  the  elongated  trough  in 
solution  space  has  disappeared  (Figure  V-3b),  and  the  region  of  low  misfit  covers  a  much 
smaller  area.  The  added  stability  from  the  regional  data  is  also  observable  in  the  inversion. 
Using  a  damping  factor  of  0.1,  a  solution  of  0  =  82°,  S  =  68°,  A  =  -89°  was  obtained  in  10 

iterations  (exact  solution:  6  =  80°,  6  =  70°,  A  -  -90°). 

In  summary,  these  test  cases  have  shown  that  the  addition  of  regional  waveforms  10  the 
set  of  observations  provides  a  valuable  constraint  for  the  inversion  of  model  parameter 
estimates  and  stabilizes  the  inversion  process. 

Advantages  of  SVD 

A  final  example  here  demonstrates  the  insight  that  SVD  provides.  This  example  is 
identical  to  that  shown  in  Figure  V-3b,  except  that  noise  was  added  to  the  radial  component 
of  the  regional  waveform  at  station  12.  With  only  one  noisy  observation,  the  final  model 
parameter  estimates  were  9  =  92°,  S  =  68°,  and  A  =  -105°.  Thus,  dip  was  unaffected  by  the 
noisy  observation,  but  strike  and  rake  deviate  by  as  much  as  15°  from  the  correct  values.  An 
examination  of  the  eigenvectors  and  eigenvalues  reveal  why  this  observation  is  having  such  a 
profound  influence  on  the  strike  and  rake.  After  eight  iterations,  only  the  noisy  observation 


is  showing  a  significant  misfit.  Due  to  the  misfit  of  the  noisy  observation,  the  error  function 
vector  has  a  large  projection  onto  the  us  data  eigenvector,  which  is  associated  with  the 
smallest  single  value.  The  weighting  factor  a3  is  therefore  large  (see  equation  2).  The  large 
factor  a3  is  weighting  the  third  model  eigenvector  v3,  which  has  components  in  the  strike  and 
rake  directions,  but  not  in  the  dip  direction.  The  parameters  strike  and  rake  are  therefore 
changing  significantly  with  each  iteration,  even  though  the  other  17  waveforms  are  being  fit. 
Thus,  the  eigenvectors  show  that  one  noisy  observation  is  imposing  a  dominant  control  on  the 
parameters  strike  and  rake,  but  is  not  influencing  dip. 

In  this  case,  the  use  of  the  maximum  likelihood  inverse  or  weighted  least-squares  is 
necessary.  Table  V-4  shows  the  inversion  results  using  different  weighting  factors  for  the 
noisy  observation.  Note  that  for  smaller  weighting  factors,  the  values  of  strike,  dip,  and  rake 
approach  the  results  of  the  noiseless  case  ( 8  =  82°,  6  -  68°,  A  =  89°)  described  in  the  previous 
section. 

Discussion 

This  inversion  method  is  not  limited  to  data  from  teleseismic  and  regional  distances 
alone.  As  long  as  the  appropriate  Green’s  functions  can  be  generated,  body  waves  for  all 
distances  can  be  inverted  for  source  parameters.  For  example,  if  the  upper  mantle  structure 
is  known,  then  waveforms  recorded  at  upper  mantle  distances  can  be  added  to  a  data  set  to 
provide  better  coverage  and  constraints.  We  have  shown  that  the  addition  of  regional 
distance  waveforms  to  the  data  set  containing  teleseismic  information,  alone,  provides  valuable 
constraints  on  the  source  parameter  estimates.  The  robustness  of  the  regional  distance 
Green’s  functions  makes  the  regional  information  valuable  for  constraining  source  parameter 
estimates,  particularly  when  the  Green’s  functions  for  the  teleseismic  portion  of  the  data  set 
are  inadequate.  The  use  of  singular  value  decomposition  for  the  inversion  is  also  particularly 
useful  since  it  allows  direct  insight  into  which  observations  are  providing  the  greatest 
influence  on  the  model  parameter  estimates. 


TABLE  V-4 


INVERSION  RESULTS  USING  MAXIMUM  LIKELIHOOD  INVERSE 
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This  method  should  also  prove  useful  for  the  inversion  of  historical  data  where 
instrument  gain  and  absolute  amplitude  are  poorly  determined  and  teleseismic  station  coverage 
is  sparse. 
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